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Abstract 

To understand the isospin effects in nearly symmetric nuclear matter we performed a complete 
pq ' systematical theoretical study within an Isospin dependent Quantum Molecular Dynamical Model 

(IQMD) and using Minimum Spanning Tree (MST) algorithm. Simulations are carried out for the 
reactions ^^^/Snso +^^*^ Sn^o and ^^"^Sn^Q +124 ^j^^q^ The collision geometry is varied from central 

^ ' to peripheral. We find that neutron rich colliding nuclei are better candidate to study the isospin 
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effects. 
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I. INTRODUCTION 



Heavy-ion collisions have been extensively studied over the last decades. The behavior of 
nuclear matter under the extreme conditions of temperature, density, angular momentum 
etc., is a very important aspect of heavy-ion physics. One of the important quantity which 



has been used extensively to study this hot and dense nuclear matter is the collective 



transverse in-plane flow 



l|. This quantity vanishes at a certain incident energy. Finite 



nuclei studies predict values for the symmetry energy at saturation of the order of 30-35 
MeV. In heavy ion collisions highly compressed matter can be formed for short time scales, 
thus the study of such a dynamical process can provide useful information on the high 
density symmetry energy. Even at low incident energies which belong to even smaller 
baryonic densities, the isospin dependence of the mean field potential was shown to yield 
same result obtained with potentials that has no isospin dependences. These results are 



in similar lines and it also indicates that even binary phenomena like fission will also be 
insensitive towards isospin dependence of the dynamics 2|]. Recently theoretical studies on 



the high density symmetry energy have been started by investigating heavy ion collisions 
of asymmetric systems 3|, \^ . Comparisons of collisions of neutron- rich to that of neutron- 



deficient systems provide a means of probing the asymmetry term experimentally 



a 



The experimental analysis of the isospin effects on fragment production has yielded several 



interesting observations: Dempsey et al. 



in their investigation of ^'^'^'^^^Xe -|-ii2,i24 g^ ^^^ 



55MeV/ nucleoli found that multiplicity of IMF's increases with the neutron excess of the 
system. A more comprehensive study was carried out by Buyukcizmeci et al. 0\ showed 
that symmetry energy of the hot fragments produced in the statistical freeze-out is very 
important for isotope distributions, but its influence is not very large on the mean fragment 
mass distributions. Symmetry energy effect on isotope distributions can survive after 



secondary de-excitation. Moreover Schmidt et al. [10| in their investigation on the analysis 
of LCP's production and isospin dependences of ^'^^Sn+^'^ Ni, ^'^^Sn+^^ Ni, ^'^^Sn+'^'^ Al at 
35MeV/nucleon and 25MeV/nucleon collisions found that isospin effects were demonstrated 
in the observables, such as the angular distribution of light particles emitted in central 



collisions at 35MeV/nucleon and LCP's emission. On the other hand Tsang et al. ll| in 
their investigation of ^^'^Sn -|-^^^ Sn, ^^^S'n-f^^^Sn systems at an incident energy of E=50 
MeV/nucleon showed the effects of isospin diffusion by investigating heavy-ion collisions 
with comparable diffusion and collision time scales. They showed that the isospin diffusion 
reflects driving forces arising from the asymmetry term of the EOS. With the passage of 
time, isospin degree of freedom in terms of symmetry energy and nucleon-nucleon cross 
section is found to affect the balance energy or energy of vanishing flow and related 



phenomenon in heavy-ion collisions [12|]. Our present study will shed light on isospin effects 
on multiplicity of fragments. We present microscopic predictions of nucleon-nucleon cross 
section in isospin asymmetric nuclear matter. In asymmetric matter, the cross section 
becomes isospin dependent. It depends on the relative proton and neutron concentrations, 
which of course also implies that the pp/nn and the np cases will in general be different 



from each other. 

This study is done within the framework of isospin-dependent quantum molecular dynamics 
model that is explained in section-II. The results are presented in section-Ill. We present 
summary in section-IV. 



II. ISOSPIN-DEPENDENT QUANTUM MOLECULAR DYNAMICS (IQMD) 
MODEL 

Theoretically many models have been developed to study the heavy ion collisions at in- 
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l|, 



termediate energies. One of them is quantum molecular dynamical model (QMD) 
which incorporates N-body correlations as well as nuclear EOS along with important quan- 
tum features like Pauli- blocking and particle production. 
In past decade, several refinements and improvements were made over the original QMD. 



The IQMD 15| model overcomes the difficulty as it not only describe the ground state prop- 
erties of individual nuclei at initial time but also their time evolution. In order to explain 
experimental results in much better way and to describe the isospin effect appropriately, the 
original version of QMD model was improved which is known as isospin-dependent quantum 
molecular dynamics (IQMD) model. 



15| model treats different 



The isospin-dependent quantum molecular dynamics (IQMD) 
charge states of nucleons, deltas and pions explicitly, as inherited from the VUU model. 
The IQMD model has been used successfully for the analysis of large number of observables 
from low to relativistic energies. The isospin degree of freedom enters into the calculations 



via both cross-sections and mean field. 

In this model,baryons are represented by Gaussian-shaped density distributions 

Mf,p,t) = J- e-(^^--W)^^ e-(-^^W)^f^. (1) 

Nucleons are initiahzed in a sphere with radius R = 1.12A^^^ fm, in accordance with the 
hquid drop modeL Each nucleon occupies a volume of h^, so that phase space is uniformly 
filled. The initial momenta are randomly chosen between and Fermi momentum(p^). The 
nucleons of target and projectile interact via two and three-body Skyrme forces and Yukawa 
potential. The isospin degree of freedom is treated explicitly by employing a symmetry 
potential and exphcit Coulomb forces between protons of colliding target and projectile. 
This helps in achieving correct distribution of protons and neutrons within nucleus. 
The hadrons propagate using Hamilton equations of motion: 

dvi d{ H ) dpi d{ H ) 



dt dpi ' dt dvi 



with 



H) = {T)+{V) 



(2) 



X fj if , p ^ t)dfdr' dpdf) . (3) 



The baryon-baryon potential V^\ in the above relation, reads as: 



V {I I ) — ^ Skyrme ^ ^Yukawa ^ ^ Coul ^ ^ sy: 



ti5{r' - f) + t25{r' - r)p 



exp{\r' -r\/p) ZiZjC^ 

+ ts -7— _, . . + 



sym 



^_i I f" + r 



(IP — r|//i) Ir' — r| 

+t,-T^TiS{r-l'-r-). (4) 

Here Zi and Zj denote the charges of i^^ and j*^ baryon, and Tg, T| are their respective 
T3 components (i.e. 1/2 for protons and -1/2 for neutrons). Meson potential consists of 

Coulomb interaction only. The parameters fi and ti, ,tQ are adjusted to the real part of 

the nucleonic optical potential. For the density dependence of nucleon optical potential, 
standard Skyrme-type parameterization is employed. The choice of equation of state (or 
compressibility) is still controversial one. Many studies advocate softer matter, whereas, 
much more believe the matter to be harder in nature. We shall use both hard (H) and soft 
(S) equations of state that have compressibilities of 380 and 200 MeV, respectively. 



The binary nucleon- nucleon collisions are included by employing the collision term of well 
known VUU-BUU equation. The binary collisions are done stochastically, in a similar way 
as are done in all transport models. During the propagation, two nucleons are supposed to 
suffer a binary collision if the distance between their centroids 



In - rj\ < J^,(Ttot = (T(^/s,type), (5) 

"type" denotes the ingoing collision partners (N-N, N-A, N-tt,..). In addition, Pauli blocking 



(of the final state) of baryons is taken into account by checking the phase space densities 
in the final states. The final phase space fractions Pi and P2 which are already occupied 
by other nucleons are determined for each of the scattering baryons. The collision is then 
blocked with probability 

Puock = 1-(1-Pl)(l-P2). (6) 

The delta decays are checked in an analogous fashion with respect to the phase space of 
the resulting nucleons. 



III. RESULTS AND DISCUSSION 

For the present study we have simulated 5oSn^'^'^+r,QSn^'^'^ and so^n^^^+so^n^^^ reactions 
by using isospin dependent quantum molecular dynamics (IQMD) model at incident energy 
600 MeV/nucleon at scaled impact parameters ( 0.0, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9). 
The collision geometry used is from central to peripheral one. The phase space generatec 



using IQMD model is analyzed by using minimum spanning tree [MST] algorithm 16 1 



16| [MSTP]. The results obtained are 



and minimum spanning tree with momentum cut 

discussed as: 

In Fig^Hshows multiplicity of free nucleons, LCP's and IMF's as a function of scaled 

impact parameters for r^oSn^'^'^ +50 Sn^^* and ^oSn^^'^ +50 Sn^^"^ . It has been observed that 

as we move from central to peripheral collisions the number of free nucleons and LCP's 

decreases because the participation zone decreases which leads to the lower number of free 



nucleons and LCP's. But in case of IMF's, curve shows a rise and fall this is because 
for central collision the overlapping of participant and spectator zone is maximum so 
we get very small number of IMF's. For semi peripheral collisions the participant and 
spectator zone decreases so the production of IMF's increases and for peripheral collisions 
very small portion of target and projectile overlap so again few number of IMF's observed 
most of the fragments goes out as heavy mass fragments (HMF's). Moreover in Fig.l 
we observe small production of free nucleons, LCP's and IMF's in ^qSu^'^'^ +50 Sn^"^^ as 
compare to 505*^^^^ +50 Sn^^*. This is because as the N/Z ratio increases the Coulomb 
repulsions increases which leads to the production of large number of products. The equation 

E{p) = E{po){-^y (7) 

Po 

gives us the theoretical conjecture of how symmetry energy varies against density. 7, tells 
us the stiffness of the symmetry energy 17|. In Fig.l the free nucleons, LCP's and IMF's for 
7 = and 7 = 0.66 both the curves clearly indicates the density dependence of symmetry 
energy. Small difference is observed in both curves in case of LCP's at a scaled impact 
parameter range from 0.0 to 0.4. Moreover when we apply MSTP cut the number of free 
nucleons increases because at low impact parameter participant zone increases and large 
number of free nucleons produced on the other hand value of LCP's and IMF's decreases 
with MSTP cut. 

Fig|2l shows the variation of multiplicity of free nucleons, LCP's and IMF's with energy 
at fixed scaled impact parameter for 5oSn^'^^ +50 Sn^^'^ and 5o5'n^°^ +50 Sn^^'^. It has been 
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FIG. 1: Multiplicity of free nucleons, LCP's and IMF's as a function of scaled impact parameter, 
observed that multiplicity of free nucleons and LCP's increases with increase in energy. On 



the other hand, one can see a rise and fall in the mu 



similar to the behaviour predicted by Aladin group 18| . Moreover number of free nucleons 



tiplicity of IMF's; this behaviour is 



and LCP's produced is very large as compare to IMF's this is because for central collision, 
collisions are violent so large number of free nucleons and LCP's produced as compare to 



IMF's. It is clear from the figure that slope of the curve is steeper in case of ^oSn +5oSn 



,124 



than 5oSn +50 Sn and this theoretical observation is in agreement with the experimen- 



tal observation of Sfienti et al. 



18| . This rise is due to the fact that in case of neutron rich 



system, heavy residues with low excitation energy will predominantly emit neutrons, a chan- 
nel that is suppressed in case of neutron-poor nuclei. We plot the curves for two different 



versions of the in-medium nucleon-nucleon cross-section term. The first one considers the 
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FIG. 2: Multiplicity of free nucleons, LCP's and IMF's with energy at fixed scaled impact parameter 

for 505'nl24 ^^^ ^^124 g^^^ ^oSn^07 +^^ 5j^124_ 



different experimental neutron-proton, neutron-neutron, and proton-proton cross- section 



and is called the isospin dependent cross section a^. r. |l5| . The second one, anoiso , considers 
the same cross section for all the isospin channels 13|, |lj]. For these two cross sections we 



observe 0.80% difference in free nucleons, 1.419% in LCP's and 1.176% in case of IMF's 
for the reaction soS'n^"'^ +50 Sn^^^ and 0.86% difference in free nucleons, 1.182% in LCP's 
and 1.88% in case of IMF's for the reaction ^qSu^'^^ +^q Sn^'^^. The difference is very small 



because the N/Z ratio is not far away from unity. But as N/Z increases the effect of different 



cross-section is clearly visible. The study on this topic is going on 19 1. 



In Figjni we have shown IMF's as a function of Zbound- The quantity Z^ound is defined 
as sum of all atomic charges Zi of all fragments with Zj > 2. Here we observe that at 
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FIG. 3: Multiplicity of IMF's as a function of Zbound- 
semi peripheral collisions multiplicity IMF shows a peak because most of the spectator 
source does not take part in collision and large number of IMF's are observed. In case of 
central collision the collisions are violent so there few number of IMF's observed and for 
peripheral collisions very small portion of target and projectile overlap so again few number 
of IMF's observed most of the fragments goes out in heavy mass fragments (HMF's). In 
this way we get a clear "rise and fall" in multifragmentation emission. It is observed that 
IMF's shows the agreement with data at low impact parameters but fails at intermediate 
impact parameters due to no acess to filters. Moreover it has been observed that for the 
isospin dependent cross-section the curve shifts towards the experimental data for both 
50^^''' +50 5ni24 and ,oSn'^' +50 Sn'^'' . 



11 



Ackno^vledgments 

This work has been supported by a grant from the university grant commission, 
Government of India [Grant No. 39-858/2010 (SR)]. 
Authors are greatful to Dr. R. K. Puri for very fruitful discussion. 



[1] A. D. Sood and R. K. Puri, Phys. Rev. C 69, 054612(2004); A. D. Sood et al, Phys. Lett. B 
594, 260(2004); A. D. Sood et al, Eur. Phys. J. A 30, 571(2006); A. D. Sood et al, Phys. 
Rev. C 73, 067602(2006); A. D. Sood et al, Phys. Rev. C 79, 064618(2009); R. Chugh and 
R. K. Puri, Phys. Rev. C 82, 014603(2010). S.Kumar, M. K. Sharma, R. K. Puri, Phys. Rev. 
C 58, 3494(1998); A. D. Sood and R. K .Puri, Phys. Rev. C 70, 034611(2004); S. Gautam, 
A. D. Sood, R. K. Puri, J. Aichehn, Phys. Rev. C 83, 014603(2011); ibid 83, 034606(2011); 
S. Goyal and R. K. Puri Nucl. Phys. A 853, 164(2011). 

[2] I. Dutt and R. K. Puri, Phys. Rev. C 81, 047601(2010); ibid 81, 044615(2010); ibid 81, 
064608(2010). 

[3] Bao-An Li, Phys. Rev. C 67, 017601(2003). 

[4] V. Greco et al, Phys. Lett. B 562, 215(2003). 

[5] W. P. Tan et al, Phys. Rev. C 64, 051901 (R)(2001). 

[6] M. B. Tsang et al, Phys. Rev. Lett. 86, 5023(2001). 

[7] H. Xu et al, Phys. Rev. Lett. 85, 716(2000). 

[8] J. F. Dempsey et al, Phys. Rev. C 54, 4(1996). 

[9] N. Buyukcizmeci et al, Eur. Phys. Journal A 25, 57(2005). 
[10] K. Schmidt et al, Acta Physica Polonica B, 41(2010). 
[11] M. B. Tsang et al, Phys. Rev. Lett. 92, 4(2004). 

[12] S. Gautam et al, J. Phys. G 37, 085102(2010).; S. Kumar, S. Kumar and R. K. Puri, Phys. 

Rev. C 81, 014601(2010); S. Kumar, S. Kumar and R. K. Puri, Phys. Rev. C 81, 014611(2010). 

[13] J. Aichelin, Phys. Rep. 202, 233(1991); R. K. Puri, et al., Nucl. Phys. A 575, 733(1994); 



12 



ibid. J. Comp. Phys. 162, 245(2000); E. Lehmann, R. K. Puri, A. Faessler, G. Batko, and S. 

W. Huang, Phys. Rev. C 51, 2113(1995); ibid. Prog. Part. Nucl. Phys. 30, 219(1993). Y. K. 

Vermani et al, J. Phys. G. Nucl. Part. Phys. 36, 0105103(2009); ibid 37, 015105(2010); ibid 

Phys. Rev. C 79, 064613(2009); ibid Nucl. Phys. A 847, 243(2010). 
[14] S. Kumar, S. Kumar and R. K. Puri, Phys. Rev. C 78, 064602(2008); ibid 57, 2744(1998); S. 

Goyal and R. K. Puri, Phys. Rev. C 83, 47601(2011). 
[15] C. Hartnack et al., Eur. Phys. J. A 1, 151(1998); S.Kumar, S. Kumar and R. K. Puri, Phys. 

Rev. C 81, 014611(2010); V. Kaur, S.Kumar and R. K. Puri, Phys. Lett. B, 697, 512(2011); 

V. Kaur, S. Kumar and R. K. Puri, Nucl. Phys. A 861 , 37(2011). 
[16] S. Kumar and R. K. Puri, Phys. Rev.C 58, 2858(1998); ibid 58, 320(1998). 
[17] H. Stocker and W. Greiner, Phys. Rep. 137, 277(1986); A. D. Sood and R. K. Puri, Phys. Rev.C 

69, 054612(2004). 
[18] C. Sfinti et. al, Acta Physica Polonica B, 37, 193(2006); ibid Phys. Rev. Lett. 102, 

152701(2006). 
[19] A. Jain, S. Kumar and R. K. Puri, Phys. Rev.C (To be submitted). 



13 



